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Field collapse, which occurs in various nonlinear systems, has attracted much attention, owing to its 
universality, complexity, and applicability. A great challenge and expectation is to achieve the controllable 
and designable collapsing pattern. Here we predict theoretically and demonstrate experimentally the novel 
collapsing behaviors of the vector optical fields in a self-focusing Kerr medium. Surprisingly, the results 
reveal that the collapse of the vector optical field is controllable and designable by engineering the 
distribution of hybrid states of polarization, and has the robust feature insensitive to the random noise. Our 
idea has its significance which it opens a new window for manipulating the optical field and the different 
kinds of field, and then facilitates to push the related researches. 

Optical field will undergo collapse at the input power above a certain critical power when an intense laser 
field propagates in a transparent self-focusing medium 15 . The high intensity as the collapse exacerbates 
will induce other nonlinear effects, such as multiphoton and saturable absorptions, to counterbalance 
self-focusing, resulting ultimately in filamentation 3 " 12 . Filamentation is an interesting and important topic due to 
its underlying physics and practical implication in many branches of physics 13 " 18 . A great challenge and expecta- 
tion is to tame the randomicity and then to realize the controllable and even designable filamentation. For this 
purpose, several methods have been proposed such as, controlling the input power and divergence angle 19 , 
shaping the field profile 20,21 , using amplitude/phase mask 22 " 25 , and introducing spatial regularization 26 . 

Polarization, as an intrinsic and fundamental vectorial nature of light, plays an indispensable role in the 
interaction of light with matter. Trisorio and Hauri 27 demonstrated the spatial control on two filaments using 
circularly polarized femtosecond laser pulses. Fibich and Han 3,28 presented a theoretical prediction that under the 
nonparaxial condition solely, the vectorial (polarization) effect of the homogeneously linearly-polarized field can 
lead to two filaments. However, polarization and its diversity of spatial distribution are not yet utilized fully. For 
the field collapse and the multiple filamentation, besides the self-focusing effect which is a requirement for the 
propagation medium, a prerequisite is that the optical field should be axial-symmetry breaking in and of itself. 
Therefore, to achieve the active control for the field collapse and the multiple filamentation, the key points are that 
the axial- symmetry breaking must be controllable and even designable, and the randomicity must be effectively 
suppressed simultaneously. A question is raised that whether the state of polarization (SoP) of light, as a degree of 
freedom, can be used to break the axial symmetry of optical field and to then availably control the field collapse 
and multiple filamentation. Vector optical fields with spatially inhomogeneous SoP distribution 29,30 provide 
indeed such an opportunity and possibility, due to the designability and diversity of SoP distribution. 

Results 

In this paper, we present an idea that the axial- symmetry breaking of optical field by engineering the azimuthally- 
variant SoP distribution controls the collapsing filamentation. We perform, to the best of our knowledge, the first 
study of the collapse of the vector optical fields in a self-focusing Kerr medium. The experimental results are in 
excellent agreement with our theoretical prediction. Some novel features are found, including the controllable, 
designable and robust features. The present study can contribute further to the study of collapsing filamentation 
in other nonlinear systems and other kinds of fields. 

An important route breaking the axial symmetry of optical field is based on the azimuthally-variant hybrid- 
polarized vector field (AV-HP-VF), which is a family of vector optical field with a "top-hat-like" intensity profile 
excluding a central singularity caused by the polarization uncertainty 30 . Under the weak focusing condition, the 
focused AV-HP-VF exhibits a uniform-intensity focal ring instead of a focal spot 30 . Its radially variant amplitude 
A(r) can be described by the generalized hypergeometric function; nevertheless, as a well approximation, A(r) can 
take the form of A(r) = a 0 r exp ( — r 1 jlr 1 ^) , where r 0 is the radius of the focal ring. Inasmuch as the weak focusing 
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does not change the SoP distribution, the focused AV-HP-VF in a 
polar coordinate system (r, <j>) can be written as 

E(r,0) = A(r) [cos S\H) + e~^ 2 sin <5| V)] , (1) 

where 3 = m§ + (j) 0 (m is the topological charge and (j) 0 is the initial 
phase), and \H) and | V) are two unit vectors indicating the horizont- 
ally and vertically polarized components, respectively. For the AV- 
HP-VF with m and (j) 0 , the local linear polarizations are located at a 
series of azimuthal positions <j> lin (n) = nnllm — (/) 0 /m, whereas local 
circular polarizations are at a series of (j) cir (n) = (2n+l)n/4m — (j) 0 /m 
(where n = 0,1, ...,4m — 1). Differently from the optical vortex 31 , the 
AV-HP-VFs carry no orbital angular momentum though they also 
have a strictly zero field at the central singularity 32 . By contrast, an 
azimuthally -variant linearly -polarized vector field (AV-LP-VF) can 
be written as follows 

E(r,^) = A(r) [cos S\H} + smS\V}}. (2) 

Compared Eq. (2) with Eq. (1), it can be found the unique difference 
that the horizontally and vertically polarized components have a 
phase difference of — n/2 for AV-HP-VF, while are in phase for 
AV-LP-VF. 

To theoretically explore the nonlinear propagation behavior of a 
vector optical field in a Kerr medium, the (2 + 1) -dimensional vector- 
version nonlinear Schrodinger (NLS) equation should be used. 
Under the slowly varying amplitude approximation, the vector-ver- 
sion NLS equation can be divided into a pair of coupled NLS equa- 
tions for two orthogonal horizontal and vertical components 

8 ~W = i v2 ^ H+ ff c ( 3 I^i 2 ^ +2 I^i 2 ^+^h)> ( 3a ) 

1% = {^y+ (3|^vlVv + 2|^ H |Vv + ^v). (3b) 
where \j/ q (q = H, V) is the nondimensional ^-component 
normalized by the total field, \j/ q (p,(j)]Q = E q (p,<j)]C) j 

y/j [|£hW; Ol 2 + \EvM\ 0\ 2 ]pdpdl p = r/r 0 , C = z/L d 
(L d = nr 0 2 /k) are the nondimensional cylindrical coordinates. 
P = 2noSoc J [\E H (r,(j)]z)\ + |£y(p,0; z)\ ] rdrd(j) is the power, P c = 
otk 2 /4nn 0 n 2 is the critical power for self-focusing, n 0 and n 2 are the 
linear and nonlinear refractive indices of the Kerr medium, a is a 
constant dependent on the initial field shape 33 , c is the speed of light 
in vacuum, and k is the wavelength, respectively. On the right- 
hand side of Eq. (3), the first term represents the contribution 
from the diffraction described by the transverse Laplacian 

= d 2 /dp 2 H — d/dp-{ — -d 2 /dd 2 , and the second term is from the 
p p 2 

Kerr nonlinearity. 

We performed a large number of simulations using Eq. (3) under P 
= 10P C in a self-focusing Kerr medium. When the spatial random 
noise is added in simulation, the random noise with a level of ± 10% 
amplitude is used. The first and second rows in Fig. 1 illustrate the 
simulation results for the collapsing behavior of the focused AV-HP- 
VF (with m = 1 and 0 O = 0) described by Eq. (1), without and with 
random noise, respectively. By comparison, the third and fourth 
rows in Fig. 1 show the simulation results for the AV-LP-VF (with 
m = 1 and <j> 0 = 0) described by Eq. (2), which is the well-known 
radially-polarized vector field (RP-VF) 29 , without and with random 
noise. Although both AV-HP-VF and RP-VF display the same focus- 
ing ring, their nonlinear propagation behaviors are quite different. 
Figure 2 shows the simulated collapsing patterns of the AV-HP-VFs 
with m = 1 (0o — 0 and n/4) and m = 2 = 0 and n/S) when the 
spatial random noise is added. Summarizing all the simulation 
results performed, we can draw following conclusions. For the AV- 
HP-VF with the given m and (j) 0 , (i) the azimuthally- variant hybrid 



C, = 0 C = 0.12 g = 0.24 £ = 0.36 




Figure 1 | Simulated nonlinear propagation behaviors of AV-HP-VF 
with m = 1 and ^ 0 = 0 in equation (2). The first and second rows 
correspond to the cases without and with random noise, respectively. For 
comparison, the RP-VFs without and with random noise are also shown in 
the third and fourth rows. The four columns (from left to right) 
correspond to four propagation distances (£ = 0, 0.12, 0.24, 0.36). The SoP 
distributions are also shown in the first column. 

SoP distribution leads to the axial- symmetry breaking; (ii) the field 
undergoes the collapse to converge into the deterministic filaments; 
(iii) the number of filaments is 4m depending solely on m; (iv) the 
filaments are always located at the azimuthal positions (j) lin (n), where 
their local SoPs are linearly polarized; (v) the collapsing pattern 
exhibits a 4m-fold rotation (C 4m ) symmetry; and (vi) the collapsing 
patterns persist among simulations, with and without different ran- 
dom noise, implying that collapsing patterns are insensitive to noise. 
Meanwhile, for the AV-LP-VFs, (i) the purely ideal AV-LP-VF 
always remains in the axially symmetric ring, indicating that the 
purely ideal AV-LP-VF cannot converge to the deterministic mul- 
tiple filaments; and (ii) the collapsing filaments produced by AV-LP- 
VF with the random noise has the uncertainty, similar with the 
collapses of optical vortices 31 and polarization vortices 34 . 

To experimentally confirm our theoretical prediction and numer- 
ical simulations, all the required vector fields with different m and (j) 0 
were created by using the experimental scheme 30 . The created AV- 
HP-VFs and AV-LP-VFs were weakly focused into a uniform- 
intensity focal ring and then incident into the CS 2 cell with a length 
of 10 mm, by an achromatic lens with a focal length of/= 300 mm. 
The intensity patterns of the field transmitted from the CS 2 cell 
are imaged on a detector (Beamview, Coherent Inc.) using an 
achromatic lens with a focal length of/ = 60 mm. The experimental 
results of the AV-LP-VFs (not shown here) reveal that the number 
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Figure 2 | Simulated collapsing patterns of AV-HP-VFs with different m 
and at a propagation distance of £ = 0.36 with random noise. 
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and positions of the induced collapsing filaments have the uncer- 
tainty among the experiments, depending on the spatial random 
fluctuation of the input field. Any one of the measured collapsing 
patterns of the m = 1 charge AV-HP-VFs, with four different 0 O 
(=0, 7r/8, te/4, 3n/S), consists of four deterministic filaments and 
exhibits a fourfold rotation (C 4 ) symmetry (Fig. 3). The four patterns 
have no difference in nature, whereas their unique difference is that 
the patterns from Fig. 3a-d are rotated counterclockwisely in 
sequence by an angle of n/S. Figure 4 shows the measured collapsing 
patterns of AV-HP-VFs (with (j) 0 = 0) for m = 1, 2, 3, 4. The m- 
charge AV-HP-VF produces 4m deterministic filaments, where m = 
1, 2, 3, 4, with the number of filaments being 4, 8, 12, 16, respectively. 
The collapsing pattern exhibits a 4m-fold rotation (C 4m ) symmetry. 
In particular, the positions of filaments for any AV-HP-VF always 
coincide with those of the local linear polarizations in the azimuthal 
dimension. The collapsing patterns of the AV-HP-VF with a given m 
and (j) 0 persist among experiments, implying that the collapsing pat- 
terns are indeed insensitive to the ransom noise. Moreover, the col- 
lapsing pattern of the AV-HP-VF is controllable and designable by 
engineering the spatial hybrid SoP structure. 

Discussion 

To understand the physics behind the collapse of AV-HP-VF, a 
cross-coupling model is presented (Fig. 5). Multiple filamentation 
requires the simultaneous spatial confinements in both the radial and 
azimuthal dimensions. The self- confinement in the radial dimension 
is due to the focal ring. However, it is of great interest and of crucial 
importance to reveal the mechanism behind the axial symmetry 
breaking that leads to the multiple filamentation. Therefore, we 
ignore the radial terms d 2 /dp 2 and d/dp in the transverse 
Laplacian V^_, while concentrate only on the azimuthal term 
d 2 /d(f) 2 in Eq. (3). The Kerr medium is divided into N equal-length 
thin segments (each thin segment has a normalized length of A£ = LI 
AT), so that the nondepletion approximation is valid within any thin 
segment. After the focused AV-LP-VF described in Eq. (2) passing 
though the 1st segment, based on Eq. (3), \j/ H and i// v can be 
expressed as follows 

^ H \ c = AC = [A+j(B + C)]cosd, (4a) 



ifr v \ c = AC = [A+j(B + C)]smd, 



(4b) 



thereby obtaining 
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Figure 4 | Measured collapsing patterns of AV-HP-VFs with 0 O = 0 for 
different m. (a) m = 1. (b) m = 2. (c) m = 3. (d) m = 4. 

where B= —Am 2 A^/A and C = aPA 3 At/P c . Based on Eqs. (4) and 
(5), we can find: (i) the horizontal and vertical components are 
always held in phase, implying that spatial SoP distribution has no 
change; and (ii) the intensity distribution still exhibits axial sym- 
metry. Therefore, a purely ideal AV-LP-VF can never produce the 
deterministic multiple filamentation, which is in excellent agreement 
with the simulation (the third row of Fig. 1). As such, any azimuthal 
position is equivalent or indistinguishable. In the presence of azi- 
muthal perturbations, axial symmetry breaking will lead to the fila- 
mentation 1 " 5 (also see the fourth row of Fig. 1). 

After the AV-HP-VF described in Eq. (1) passing through the 1st 
segment, based on Eq. (3), \j/ H and xj/y can be expressed as follows 

^h\ C =ac = i A +A B + ( C / 3 ) C 1 + 2 cos2 S )}} cos S > ( 6a ) 

^ v \^ A ^=j{A+j[B+{c/?>)(i + 2sm 2 5)}}sm5, ( 6b ) 

thereby obtaining 

(|M 2 + |^| 2 )| c= , c = A'+B'cos(4<S) 



where 



A'=A 2 + (B — C) - (B + 5C/6) (C/3), 
B'=-(B-5C/6)(C/3). 



(7a) 

(7b) 
(7c) 




Figure 3 | Measured collapsing patterns of AV-HP-VFs with m = 1 for 
different 0 O . (a) 0o = 0. (b) 0 O = nIS. (c) 0 O = n 14. (d) 0 O = 3tc /8. 



0.2 0.4 0.6 0.8 
Azimuthal position <j) (n) 

Figure 5 | Cross-coupling model for the collapse of AV-HP-VF. (a) The 

intensity and SoP distributions in the input plane, (b) Those behind the 
first segment, (c) Those behind the rctfi segment, (d) Those in the output 
plane. Arc arrows indicates the energy exchange directions, (e) The 
simulated evolution of S 3 of azimuthal hybrid SoPs for the AV-HP-VF with 
m = 1 and 0 O = 0 (only the azimuthal range from (j) = 0 to n is shown). 
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Figure 6 | Azimuthal self-focusing model for the collapse of AV-HP-VF. (a) Input focused AV-HP-VF. (b) Induced azimuthal self-focusing lenses, 
(c) Axial- symmetry broken AV-HP-VF. (d) Induced azimuthal self-focusing lenses with the shorter "focal length", (e) Collapsing filaments. 



Equation 7 suggests that the nonlinear cross -coupling results in the 
energy redistribution of AV-HP-VF in the azimuthal dimension (or 
the axial symmetry breaking). It should be pointed out that the 
intensity distribution still holds the axial symmetry in the absence 
of nonlinearity (n 2 = 0— >C = 0->-B' = 0). As an example, the AV-HP- 
VF with m = 1 and <j> 0 = 0 is considered (Figs. 5a-d). In the input 
plane (£ = 0), the focused AV-HP-VF has the azimuthally-invariant 
intensity distribution, while the azimuthally- variant hybrid SoP dis- 
tribution (Fig. 5a). After the 1st segment, the cross-coupling results 
in the energy exchange in the azimuthal dimension toward the loca- 
tions at (j) lin (n) from those at ^ dr (n) (Fig. 5b). As the propagation 
distance increases, more energy is transferred toward the locations at 
(j) lin (n) after the jth segment (Fig. 5c). In the output plane (Fig. 5d), the 
cross-coupling leads to the energy almost converging into the azi- 
muthal locations at (j) lin (n). The axial symmetry of the focused AV- 
HP-VF is broken ultimately to form four deterministic filaments. 
From Eq. 6, the cross-coupling acquires an additional phase, which 
makes the phase difference between the horizontal and vertical com- 
ponents to not always preserve n/2, resulting in the SoP evolution in 
contrast to AV-LP-VF. However, at two series of azimuthal locations, 
(j) lin {n) and (j) cir {ri), the SoPs have no change. The simulated evolution 
of the azimuthal-variant hybrid SoPs, by the Stokes parameter S 3 
characterizing the ellipticity 30 , during the nonlinear propagation of 
the AV-HP-VF with m = 1 and (j) 0 = 0 confirmed the above predic- 
tion (Fig. 5e). As propagation distance increases, the local ellipticity 
has a tendency to decrease (i.e, the local polarization ellipse will 
become flatter). Except for at the special azimuthal locations of 
(j) lin (n) and (j) cir (n), the SoPs have indeed no change, like the above 
discussion from the cross-coupling model. 

To understand the collapsing filaments, the azimuthal self-focus- 
ing model is also very beneficial (Fig. 6). The self-focusing phenom- 
enon originates from the spatially variant refractive index change. In 
an isotropic Kerr medium under the same intensity, the induced 
refractive index change depends on the SoP of optical field (An lin 
> An eU > An cir ), where An lin , An el \ and An cir are the refractive index 
changes induced by the linearly, elliptically, and circularly polarized 
light, respectively 35 . As a result, the refractive index change induced 
by the AV-HP-VF is azimuthally variant, leading to the axial-sym- 
metry breaking. The local linear and circular polarizations are 
located at <j> Un (n) and (j) cir (n) (where n = 0, 1, 2, 3) for the AV-HP- 
VF with m = 1 and <j> 0 = 0 (Fig. 6a). Therefore, the induced four 
azimuthal self-focusing lenses have their centers located at <j) lin {ri) and 
their edges located at (j) cir (n) (Fig. 6b). The input field undergoes the 
azimuthal multiple-nucleation self-focusing, making the energy con- 
verge toward the four locations at (j) lin (n) (Fig. 6c). As the propagation 
distance increases, the "focal lengths" of the induced azimuthal self- 
focusing lenses become shorter (Fig. 6d), and then the energy con- 
verges ultimately into the four deterministic filaments (Fig. 6e). In 
particular, An lin is greater than 50% An ci \ which is the reason why 
AV-HP-VF can suppress the random noise (the robust feature) dur- 
ing the deterministic collapsing filamentation. 

We have presented the theoretical prediction and conducted the 
experimental evidence on the novel collapsing behaviors of AV-HP- 
VF in the self-focusing Kerr medium. The azimuthally-variant 
hybrid SoP structure is crucial in the axial symmetry breaking. The 
collapsing filamentation has the controllable and robust features. 



The m-charge AV-HP-VF produces 4m filaments with 4m-fold rota- 
tion (C 4m ) symmetry. The number of multiple filaments is a function 
solely of m, and their locations are determined by both m and <j> 0 . Our 
results suggest that engineering the hybrid SoP structure can tame 
the randomicity and then realize the designable collapsing filamen- 
tation. This study only concentrated on the collapsing behavior; thus, 
the higher-order nonlinearity, which is important in the postcollapse 
process, has been neglected. Of course, we also simulated the mul- 
tiple filamentation when the higher- order Kerr effect is considered, 
the results indicate a novel feature that the relative phase difference 
among the multiple filaments is stable and insensitivity to the input 
fluctuations of the pulse energy, which is quite different from the 
reported fact that the relative phase difference is instable and sens- 
itivity to the small fluctuations of the pulse energy, i.e., loss of phase 
of collapse 36 . Our idea may offer an alternative route to produce the 
controllable and robust multiple filamentation in other nonlinear 
systems, thereby facilitating the development of additional surpris- 
ing applications. This work has the great significance that it not only 
opens a new window for manipulating the optical field and facilitates 
to push the related researches, but also has a certain reference value 
for manipulating the different kinds of field. 

Methods 

Generation of femtosecond vector fields. We follow the same method as Ref. 30 for 
creating the vector fields required in this study. The used light source is a Ti: sapphire 
regenerative amplifier femtosecond laser system (Coherent Inc.), which provides a 
fundamental Gaussian mode with a central wavelength of 800 nm, a pulse duration of 
35 fs, and a repetition rate of 1 kHz. An achromatic 1/2 wave plate and a broadband 
polarized beam splitter were used to control incident laser fluence on the sample 
surface. Another achromatic 1/2 wave plate was used to change the polarization 
direction of the incident femtosecond laser into the vector field generation unit, where 
all the elements are achromatic to suppress pulse broadening as much as possible. The 
created AV-HP-VFs and AV-HP-VFs have a "top-hat-like" spatial profile, excluding 
a central singularity caused by the SoP uncertainty, with a radius of 1.5 mm, a pulse 
energy of 6.0 uj, and a pulse duration of —65 fs. 
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